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Strong Shock Waves and Nonequilibrium Response in a One-dimensional Gas: 

a Boltzmann Equation Approach 
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We investigate the nonequilibrium behavior of a one-dimensional binary fluid on the basis of 
Boltzmann equation, using an infinitely strong shock wave as probe. Density, velocity and temper- 
ature profiles are obtained as a function of the mixture mass ratio a. We show that temperature 
overshoots near the shock layer, and that heavy particles are denser, slower and cooler than light 
particles in the strong nonequilibrium region around the shock. The shock width which char- 

acterizes the size of this region, decreases as ~ /i 1 ^ 3 for fi —* 0. In this limit, two very different 
length scales control the fluid structure, with heavy particles equilibrating much faster than light 
ones. Hydrodynamic fields relax exponentially toward equilibrium, cj>(x) ~ exp[— x/X\. The scale 
separation is also apparent here, with two typical scales, Ai and A2, such that Ai ~ /1 1 / 2 as /1 —* 0, 
while A2, which is the slow scale controlling the fluid's asymptotic relaxation, increases to a con- 
stant value in this limit. These results are discussed at the light of recent numerical studies on the 
nonequilibrium behavior of similar f d binary fluids. 



I. INTRODUCTION 

When a piston moves at constant velocity into a fluid, 
it generates a shock wave.[l| This is a strong perturba- 
tion that drives the system far from equilibrium. Study- 
ing the structural properties of the shock wave, and 
how equilibrium is restored behind it, we may extract 
valuable information on the fluid's transport proper- 
ties, the damping of nonequilibrium fluctuations, etc. 
That is, we may use the shock wave as a probe to bet- 
ter understand the nonequilibrium response of the fluid. 
This is particularly appealing for one-dimensional (Id) 
systems. Their nonequilibrium behavior has attracted 
much attention during the last years. 0-^3 This fol- 
lows from two main reasons. On one hand, the di- 
mensional constrain characteristic of Id systems plays 
an essential role in many real structures, ranging from 
carbon nanotubesQ to anisotropic crystals 0, magnetic 
systemsjjl solid copolymers 0, semiconductor wires 
zeolites 9], Bosc-Einstein condensates [lOl, colloids in nar- 
row channels DNA, nanofluids, etc.|^J On the other 
hand, the simplicity and versatility of Id models allow 
one to tackle fundamental questions in nonequilibrium 
statistical mechanics, as those related to irreversibility, 
normal transport, equipartition, local thermodynamic 
equilibrium, etc.|2j 

In this way, it has been found that the low dimensional- 
ity can give rise to anomalous transport properties in Id 
fluids. To be more specific, the single-file (SF) constraint 
characteristic of Id fluids introduces strong correlations 
between neighboring particles, which asymptotically sup- 
press mass transport 111 and enhance heat conduction 
El El El El 111 Ellli EpL In particular, it is cur- 
rently believed (see however |l2tll^|) that Id fluids with 



momentum-conserving interactions and non-zero total 
pressure exhibit a thermal conductivity k that diverges as 
k ~ IP when the system size L^ooElElElEllEi 
However, there is no complete agreement yet on the expo- 
nent 7 > 0: there exist strong numerical and theoretical 
evidences pointing out that 7=1/3 0,0, although 
mode-coupling theories and Boltzmann equation analy- 
sis predict 7 = 2/5 lla ll 71 . and other numerical results 
are closer to 7 = 1/4 Conservation of momentum 

seems to be crucial; as soon as this symmetry is bro- 
ken, normal heat transport is recovered. 19| Moreover, 
momentum-conserving Id fluids typically exhibit other 
anomalies when perturbed away from equilibrium, as 
for instance dynamic non-equipartition of energy, emer- 
gence of multiple relaxation scales, energy localization, 
etc. El El El E El In general, the behavior of Id fluids 
far from equilibrium still poses many intriguing questions 
which require further investigation. 

In this paper we probe the nonequilibrium response 
of a Id model gas. In particular, we address the prob- 
lem of an infinitely strong shock wave propagating into 
a Id binary fluid, within the context of Boltzmann 
equation. jH HI |il| The shock wave characterizes the 
transition between two different asymptotic equilibrium 
states of the fluid. For strong shock waves (to be spec- 
ified latter), this transition happens within molecular 
scales, so kinetic (Boltzmann) equations must be used. 
In fact, one may think of the shock wave problem as 
the simplest case study dominated by the nonlinearity 
of Boltzmann equation. [23 As we will see below, Boltz- 
mann equation provides a powerful tool to investigate 
the observed structural and relaxational anomalies in Id 
fluids. Kinetic theory methods have been applied with 
success to study strong shock waves in high-dimensional 
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simple fluids. [2J, |25|, |26| However, in one dimension, a 
simple (i.e. mono-component) fluid constitutes a patho- 
logical limit. [23 This follows from the fact that parti- 
cles with equal masses in Id interchange their velocities 
upon collision, so a Id mono-component fluid can be re- 
garded to a large extent as an ideal gas of non-interacting 
particles. [2^ We may restore ergodicity and relaxation in 
velocity space by introducing different masses. The sim- 
plest case is that of a binary mixture. In particular, we 
study in this paper a binary Id fluid composed by two 
species of hard-point particles, with masses mi < m 2 and 
equal concentrations. Energy transport in this model has 
already been extensively studied via molecular dynamics 
simulations. El El E H3 

The structure of the paper is as follows. In section II, 
we write down the relevant Boltzmann equation for our 
Id binary fluid, and extend a method by H. Gradf^Efij] 
to study the structure of an infinitely strong shock wave 
propagating through the fluid. The results of this ap- 
proach are described in section III. In particular, we 
study there the shock hydrodynamic profiles, paying spe- 
cial attention to the shock width and the relaxation to- 
ward equilibrium, as a function of the mass ratio. Finally, 
in section IV we discuss our results. 



II. THE SHOCK- WAVE PROBLEM 

The Boltzmann equation for a Id binary fluid in a ref- 
erence frame moving with the shock wave reads. |2ll E2I 



= Qi(/i,/ 2 )-/iK/2), 
= Q 2 (/ 2 ,/i)-/ 2 K/i)- 



(i) 

(2) 



Here fi(x,v) is the probability density for finding a par- 
ticle of type i = 1,2 (i.e. mass m^) at position x with 
velocity v, and Qi and v represent the gain term and 
collision frequency in the collision operator, respectively. 
They are defined as, 

/+00 
dw \v-w\fi(x, v^jfj (x, w'i), (3) 
-OO 



+00 



dw \v — w\fi(x, w). 



(4) 



Index j refers here to particle species other than i, and 
(v'^Wj) are pre-collisional velocities that after collision 
give rise to velocities (v,w) for the pair namely, 



v[{v,w) 
w'i{v,w) 



(rrij — frij)v + 2m jiu 

mi + rrij 
(rrij — mi)w + IrriiV 

m, + rrij 



Notice that eqs. 0-© only include the effect of cross- 
collisions between the different species, and no self- 
collisions between like particles. This is reflecting the 
one-dimensional character of our model fluid, and it is a 
main difference with Boltzmann equation for mixtures in 
higher dimensions. 

The structure of the shock wave can be deduced from 
eqs. ©-10 subject to the boundary conditions, 



fi(x,v) -> Gi ±(v) 



where 



Gi,±(v) = n± 



2nT± 



exp 



±00, 



rriiiv - u±) 2 



2T+ 



(6) 



(7) 



(5) 



Here n± , T± and u± are respectively the number density, 
temperature and macroscopic velocity at x — ► ±oo.|3J| 
One usually chooses u + > 0, which implies (see below) 
U- > 0. This amounts to represent a flow of the fluid 
mixture from — 00 (upstream or pre-shock region) to +00 
(downstream or after-shock region) . 

We are particularly interested in the structure of an in- 
finitely strong shock wave. In this way we can probe the 
response of the Id binary mixture in the strong nonequi- 
librium regime, far from the linear response region. [22] 
The strength of a shock may be defined in several ways. 
One of them is based on the ratio of downstream to up- 
stream pressures, r = p + /p-, where p± = n±T±. Weak 
shocks have a ratio r close to unity. In this case, the shock 
thickness is large as compared to the mean free path, and 
continuum approximations of the type of Navier-Stokes 
(or Euler, Burnett, etc.) may be used to characterize the 
shock structure. However, as r grows the shock thick- 
ness becomes comparable to the mean free path, making 
inappropriate the application of hydrodynamic approxi- 
mations. In this case kinetic equations, as those in © 
|2"|. must be used instead. We pay attention in what 
follows to the limit r — > 00, or equivalently T_ — > 0. 

The parameters entering the two limiting Maxwellians 
are related due to conservation of particle number, to- 
tal momentum and total energy. Such conservation laws 
can be expressed as the constancy of the corresponding 
fluxes, which yields, 

n+u+ = n_u_, (8) 
n+ [(mi + m2)u\ + 2T + ] = n - [i m i + m 2)"_] , 
n + [(mi + m 2 )u 3 + + 6u + T + ] = n_ [(mi + m 2 )u 3 _] . 

These are the Rankine-Hugoniot conditions, which relate 
the downstream and upstream values of the flow fields in 
the mixture. [1( Solving this system for the downstream 
asymptotic flow fields one finds. |3~0] 

u _ \ 
n + = 2n_, u + = — , T + = -(m 1 +m 2 )u 2 _. (9) 
z o 

Therefore, the downstream asymptotic density (velocity) 
is twice (half) the upstream one for r — > 00. It must 
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be noticed that these results come from the conservation 
laws characterizing the two-component fluid, and have 
nothing to do with the hydrodynamic behavior of the 
mixture. In this way it is not strange that even in non- 
hydrodynamic cases, as the equal mass version of our 
one-dimensional gas,[2£| the results © hold. 

The case of a infinitely strong shock has been studied 
in detail for mono-component gases in dimensions larger 
than one EiL Efij | Extending a successful hypothesis by 
H. Grad,[24j we suggest here that the problem of a shock 
wave in our Id mixture is in general well-defined in the 
limit r — > oo, and that in this case the distribution func- 
tions fi can be decomposed into a singular part, fi, pro- 
portional to a Dirac ^-function on the upstream velocity 
U-, and a regular part, fi, smoother than the former. 
This is because for r — > oo, particles in the upstream 
region exhibit a well-defined fixed thermal velocity when 
measured on the downstream velocity scale. Hence, 

fi = fi + fi, with fi(x, v) = hi(x)S(v - uJ). (10) 

The singular particle densities hi(x) control how impor- 
tant the singular component is at a given coordinate x. 
In particular, it is evident that fn(x) — > as x — > — oo, 
and hi(x) — > as x — » +oo. Introducing the above de- 
composition <|10|) into eqs. we may split the orig- 
inal Boltzmann equations in two different sets, attending 
to the singular or regular character of the terms, 



N'(x,n) 



+ Mi) = Qiifl fi) ~ Mi), (11) 



} dfi 
dx 



Qi(kfi)-Qi{fi,fi) + fiv(fi) = 

= Qi(fi,fi)-Mi), (12) 



with i = 1,2. Eq. (|12|l for the regular component fi looks 
more complicated than the original Boltzmann equations 
CJ (J2J. However, as advantage, it lacks the singular be- 
havior characterizing for r — > oo the whole distribution 
fi. This allows us to perform a simple local Maxwellian 
approximation for fi that yields meaningful results. In 
particular, we now assume, 



fi(x,v) = hi(x) 



2TtT t (x) 



cxp 



rrii(v — u(x)Y 



2Ti(x) 



(13) 



Here hi{x), u{x) and Ti(x) are the number density, ve- 
locity and temperature associated to fi. Now hi(x) — > 
as x — y — oo and hi(x) — > n+ as x — > +oo. We have also 
assumed at this point that the regular flow velocity u(x) 
does not depend on particle species. This assumption 
seems natural given the one-dimensional character of the 
system, which forces neighboring particles to move co- 
herently on average. Therefore the problem of the shock 
wave structure reduces to computing the spatial depen- 
dence of 7 different flow fields, namely hi{x), hi(x), Tiix) 
and u(x), i = 1,2. 




FIG. 1: Total density profiles (normalized) as a function of 
mass ratio /i. Notice the logarithmic scale in the /i-axis. The 
inset shows the same profiles in a 2d plot for better compari- 
son. 




FIG. 2: Excess density profiles (normalized). 



It is important to notice that the Maxwellian approx- 
imation HI 3(1 is not a local thermodynamic equilibrium 
(LTE) hypothesis, since it only affects the regular com- 
ponent fi(x,v) of the distribution function. The singu- 
lar, delta-like component fi(x,v) accounts for the most 
important nonequilibrium effects. However, it is natural 
to question the validity of 1(13(1 for the present strong 
nonequilibrium problem. This approximation is just 
the zeroth-order term in a formal expansion of fi(x,v) 
around local Maxwellian behavior, in the spirit of Sonine 



polynomials expansions, Gram-Char lier and Edgeworth 
series, etc. l34l | One would expect that higher-order 

terms in this expansion could be relevant near the shock 
layer. In fact, such slight deviations have been observed 
in mono-component gases in higher dimensions. |26l| How- 
ever, these small corrections have negligible effects on the 
shock hydrodynamic profiles. |2q and therefore are not 
important for our discussion here. 

We may use l(l(J(l and ((13(1 to compute the local hy- 
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FIG. 3: Flow velocity profiles (normalized). 



FIG. 4: Excess flow velocity profiles (normalized). 



drodynamic currents. Let J^„(x), Ji >v [x) and Ji, e [x) be 
the particle, momentum and energy fluxes of species i at 
position x, respectively. They can be written as, 



J i,v (^) 



dv vfi{x,v) — hiU + fiiU^, (14) 



dv miV 2 fi(x,v) 
" u 2 + hiU 2 _} + fiiTi, (15) 
dv rriiV 3 fi(x,v) 
mi(niU 3 + hiii 3 _) + 3uniT i} (16) 



mi{riiU 



(A trivial factor i has been omitted in the definition of 
the energy flux for convenience). Given the constancy 
of the total fluxes along the system, we now write the 
Rankine-Hugoniot conditions as, 

fi\u + fi\U- — 7i_u_, (17) 

hiU + fi2U- = ri-U-, (18) 
mi{h\u 2 + fiiu 2 J) + m2(fi2U 2 + fi2U 2 _) + 

+Pi +P2 = n_(mi + m2)u_, (19) 
mi(hiu 3 + fiiv?_) + m2{fi2U 3 + fi2U 3 _) + 

+3u(pi +P2) = ra_(mi + m 2 )u 3 _, (20) 

where pi = h{Fi is the pressure of the regular component 
Defining M = m,\h\ + rri2n2, M = ra\h\ + rri2n2, 
M_ = (mi + rri2)n^ 1 and P = p\ + P2, we obtain, 

Mm = Au_, Mil 2 + P = Au 2 ^, Mu 3 + 3uP = Au 3 _, 

(21) 

where A = M_ — M. Using the first relation into the 
second equation one has P = Au_(u_ — it), and insert- 
ing this expression into the third equation yields (after 
division by Ait-) 2u — 3u-u + u 2 _ = 0. This equation has 
solution u = U-/2 (the second root, u — is trivial 
and can be neglected as incompatible with the boundary 



condition at +00). Therefore the regular flow velocity is 
constant along the system, and equal to its asymptotic 
value at +00, u+. Using this in eqs. (|17|I - (|18JI one has, 



hi{x) = 2[n_ 
u{x) = \u- 



i(x)], 1 = 1,2 (22) 
(23) 



Combination of these results with P = Au_(m_ — u) 
yields a relation for the regular pressures, namely, 



Pi +P2 



— (minx + m2n 2 ) 



(24) 



These manipulations are intended to express all re- 
maining flow fields in terms of the singular particle den- 
sities, hi{x). Equations for these singular densities can 
be derived integrating with respect to v eq. (|TTT) for the 
singular component £■ . |35l| This yields, 



w_ — = 2aj(x)ni{x)\rij{x) — n_ 

dx 



(25) 



where, 



2T l (x) 2 




2T l {x) 



-exp 



rriiuz 



8T,(x) 



(26) 



and where we have used results ()22|l and (|23|) . To com- 
pletely specify the flow fields, we need another expression 
relating the regular pressures pi to the other fields. This 
relation can be derived studying the momentum and en- 
ergy transfer from species i to species j. Multiplying 
eq. by rriiV (momentum transfer) or rriiV 2 (energy 
transfer) and integrating over v, we arrive to differential 
equations for Ji, v (x) and Ji >e {x) of the form, 



dJ l:V 
dx 

dJ t , e 
da; 



IL», 
(a), 



(27) 
(28) 
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where the fluxes are defined in eqs. (| 1 5 1> — (| 1 €>|) and, 

/ + OC 
d«m i »[e i (/ i ,/ i )-/ i i/(/j)], (29) 
-00 
r+oc 

Ei(x) = / d^ m^ 2 [Qi(fi, fj) - f iV (Jj)] . (30) 



Since the total momentum flux, Ji_ v {x) + J2,v( x )i an d 
total energy flux, J l e (a;) + J2, e {x), are constants along 
the line (Rankine-Hugoniot conditions), the above differ- 
ential equations, l|27l) - (|28[) . express the transfer of mo- 
mentum and energy from species i to species j at a given 
point x. Using the expressions for Ji jV (x) and Ji^ e (x) in 
(I15|l Ill6() . it is easy to show that 3uJi tV (x) — Ji : e{x) — 



^rriin^u^ — constant, and therefore, 



3uIk(x)-Ei(x) = 0. 



(31) 



This last equation will give us the desired extra rela- 
tion for the regular pressures. In general, the integrals 
in (|29|) - (|30[) cannot be performed due to the velocity- 
dependent collision kernel \v — w\ appearing in the defini- 
tion of Qi and v, see (01 and (@J. In order to continue our 
derivation, we now approximate this kernel by a generic 
velocity-independent kernel u{x), much in the spirit of 
Maxwell molecules. 29J 

Using this assumption we can solve the above integrals, 
obtaining, 



Tli(x) 
Ei(x) 



2m,m,- 



n-U-a(x)(N i -N j ), (32) 

m i m j f 2 r/ 

7 ■ -^<7{x)< u_ mi + 2mj)n i n ] 

(m l + rrij) 2 I 

(rrij + 2m l )n i h J ] - 4[Njpi - Nipj] j, (33) 



where we have defined the species total number density, 
Ni(x) = hi(x) + fii(x) = 2n_ — fii{x), i = 1,2. Within 
this approximation, eq. 13111 reduces to, 



Njpi - Nipj 



-7r( m j 



m i ){n i n : j + hifij), (34) 



which is the second relation between the regular pressures 
we were looking for. Together with l|24|l . this yields, 



u 2 _ [(rrij — m,i)(hihj + hihj) + 2MNi 



(35) 



where we have used the previous definition M = rmhi + 
nijfij. This result, based on the Maxwell velocity kernel 
approximation, gives predictions for the regular pressures 
which are very close Vx to the results obtained from the 
solution of eq. (|31(l based on the numerical evaluation of 
integrals (f2T)|l - lf3T)| . thus validating our approximation. 
The reason behind this good agreement is that Eq. Ij34(l 
does not depend on the collision kernel cr(x) chosen for 
the Maxwell-like approximation: the same equation is 



T(x,n) 




FIG. 5: Temperature profiles (normalized). 




FIG. 6: Excess temperature profiles (normalized). 



found for any kernel ct{x) one may think of, and in par- 
ticular for the (unknown) optimal kernel which better 
approximates the real velocity-dependent kernel. 

Recalling that hi = 2(n_ — hi) and u — \u-, eqs. I|35f) 
express the regular pressures pi — hiTi, and hence the 
regular temperatures, in terms of the singular densities 
hi. Therefore the problem of the shock wave reduces to 
solving the coupled, strongly non-linear differential equa- 
tions for hi(x), i = 1,2, since all other non-constant 
flow fields, hi and Ti, with i = 1,2, can be written in 
terms of fi, . 



III. RESULTS 

The strong non-linearities present in eqs. (|25|) prevent 
any full analytical solution to the shock wave problem, 
so numerical integration must be used instead. This pro- 
cedure yields the singular densities hi{x), from which 
all hydrodynamic flow fields follow. We are particu- 
larly interested here in the total number density N(x), 
flow velocity u(x), and temperature T(x), and their par- 
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ticularization to each species (Ni(x), Ui(x) and Tj(x), 
i = 1,2, respectively). The total number density is 
defined as N(x) = Ni(x) + N 2 (x), where N^x) = 
fii(x) + fii(x) has been defined above. The total flow 
velocity and temperature fields are defined respectively 
asu(x) = [J2i m i N i( x )l— J2i miNi(x)ui(x) and T(x) = 
N(x)~ l 'E i N i (x)T i (x),m where > 



-Total 



Ti 



N: 



N7 



+00 



dv vfi(x,v) 



N, 



+ OC 



dv n%i [v - u(x)] 2 fi(x, v) 



^ Uifi + mi[fii{\u^ - u) + fn{u- - u) ] 
~ Nt ' 

For a generic hydrodynamic flow field 4>(x), we now 
define an associated normalized field as <fr'(x) = {(j)(x) — 
4>-)/(4>+ — </>_), where <p± — <j>(x — > ±00). This allows 
us to compare shock profiles for different mass ratios /1 = 
mi/m 2 € (0,1). In addition, we define the normalized 
excess flow fields as 5<j>'(x) = (f> 2 (x) — 4>'\(x) to study 
the species dependence of the profiles. Our results for 
N'(x), u'(x) and T'(x) are given in Figs. 1, 3 and 5, 
respectively, while the excess fields SN'(x), 5u'(x) and 
ST'(x) are depicted in Figs. 2, 4 and 6. It should be 
noted that we set the coordinate origin x — so that 
N'(x = 0) = \. 

A general observation is that all the profiles get 
steeper, and the excess fields get more localized and 
peaked around x — 0, as n decreases. However, there 
are fundamental differences for the different hydrody- 
namic quantities. The density flow field, N'(x), con- 
verges toward a limiting shape as fi —> 0, characterized 
by a sudden increase of density from N'(x — > 0_) = to 
N'(x — > 0+) = i, followed by a much slower relaxation 
toward the asymptotic value N'(x — > +00) = 1, see Fig. 
1. On the other hand, the flow velocity profile u(x) con- 
verges toward a step-like function localized at x = 0, 
Fig. 3, while the temperature field, Fig. 5, exhibits an 
overshoot which sharpens and increases as /x — > 0. In 
addition, all profiles are markedly asymmetrical. Fig. 7 
exhibits all fields for several mass ratios. 

The excess profiles contain the information about the 
relative local distribution of the hydrodynamic fields be- 
tween the species. In particular, SN'(x) shows that the 
number density of heavy particles (type 2) around the 
shock layer is larger than that of light particles (type 1), 
see Fig. 2. This excess concentration of heavy particles 
around the shock is asymmetric, with most excess heavy 
particles lagging behind the shock, and it increases as 
fi — > 0. On the other hand, the flow velocity u\(x) of 
light particles around the shock is larger than the ve- 
locity of heavy particles, see Fig. 4, the difference in- 
creasing with decreasing fi (notice here that Su'(x) > 
involves u\(x) > 1*2(2;)). Finally, the excess tempera- 
ture field exhibits an interesting behavior, see Figs. 6 




FIG. 7: (Color online) Density, velocity and temperature pro- 
files for each species and the mixture. From left to right 
and top to bottom, fj, = 0.000125,0.001,0.016,0.128. Veloc- 
ity (temperature) profiles have been shifted one (two) units 
in the vertical axis for visual convenience. 



and 7. First, 5T'(x) is mostly negative near the shock, 
so heavy particles are cooler than light particles around 
the shock, despite the former are denser than the latter 
around x — 0. More in detail, 5T'(x) exhibits an oscil- 
lation around x = 0, meaning that heavy particles are 
much colder than light particles behind the shock layer 
(the cooler the smaller /i is), while the situation is re- 
versed in front of the shock, where heavy particles are 
now more energetic than light ones, see Fig. 6. 

To better understand some of these observations, we 
can analyze the asymptotic behavior of eqs. (|25() . First, 
we study how the downstream binary fluid goes toward 
equilibrium by letting x 3> 0. In this case, the singular 
densities are very small, > 0) C 1, so we can lin- 
earize eqs. ()25fl. obtaining a simple decoupled system of 
differential equations, 



ii- — — — = — 2aj(+oo)n-ni(x), 



(36) 



with i = 1,2, and where, 



a, (+00) = — 



mi + mi' 



mi + m 2 



-exp( 



mi +m 2 ' 



(37) 



is the limit of the amplitudes H26f> as x — > +00. There- 
fore, we find fii{x) ~ exp[— x/Xj] in this limit, with 
Xj = M_/[2n_aj(+oo)], so the downstream binary fluid 
reaches equilibrium exponentially fast, as expected from 
the uncorrelated nature of Boltzmann equation. [2ll EM 
l23| This equilibration process is characterized however by 
two different typical scales: Ai controls how heavy parti- 
cles (type 2) relax asymptotically, while X 2 controls the 
behavior of light particles (type 1). The /x-dependence 
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<, i 




FIG. 8: (Color online) Main: Mass ratio dependence of the 
typical scales Ai controlling the exponential relaxation of the 
downstream binary fluid toward equilibrium. Notice that Ai 
scales as /x 1 ^ 2 for small fi. The solid dot signals the singular 
limit jit = 1, for which both typical scales formally diverge. 
Inset: Mass ratio dependence of the shock width uj in log-log 
scale. The width decreases as w ~ ft 1 / 3 for /j, — > 0. 



of both typical scales is depicted in Fig. 8. Here we 
see that A 2 > Ai V/i 6 (0, 1). In fact, A 2 converges to 
a constant as /i — > 0, while Ai goes asymptotically to 
as /i 1 / 2 , see eq. <|37|) . Therefore, a separation of scales 
emerges for very small mass ratios /i, with both fast and 
slow evolution scales controlling the fluid relaxation. All 
the macroscopic hydrodynamic fields <fi(x) relax asymp- 
totically as exp[— X/X2], following the slowest relaxation 
scale associated to light particles, except for the number 
density for heavy particles, which relax much faster, as 
exp[— x/Xi]. 

In the opposite limit, fi — ► 1, a singularity emerges. 
In fact, Boltzmann equation does not yield any useful 
information for ji = 1 in one dimension (Id): the gain 
and loss terms in the collision operator are equal in this 
case, so the collision term in Boltzmann equation is zero 
(see right hand side of Eqs. (!}-©). That is, collisions 
do not change the velocity distributions fi(x,v) in this 
limit, as expected from the fact that equal-mass particles 
exchange their velocities upon collision in Id. Therefore 
there is no relaxation of hydrodynamic fields for /i = 1 , 
and the relaxation scales Ai are not defined. Alterna- 
tively, this can be interpreted as a divergence of both 
A,;. However, this is a singular behavior which emerges 
only for fx = 1. As soon as fj, < 1, normal behavior is 
recovered. In particular, for /i slightly smaller than 1, 
one expects both typical scales to be almost equal, as 
observed in Fig. 8. 

The separation of scales described above for /i — » 
does not only emerge for x — > +00, but characterizes the 
profiles \fx. In fact, the amplitudes cti(x), see eqs. (|26|) . 
are such that a\(x) 3> 0:2(2;) as /i — > 0. Hence, in this 
limit the singular density 77,2(2;) changes at a much faster 



rate than 711(2;) does, see eqs. (|2"5)l . so we can consider 
h\{x) as a constant in the scale in which hiix) evolves. 
This effectively decouples again the differential equations 
Q25JI (though they are still highly non-linear), giving rise 
to a multi-scale relaxation process: (i) In a first step, 
heavy particles very quickly relax toward their equilib- 
rium downstream asymptotic state while light particles 
remain close to the upstream state, (ii) This is followed 
by a relatively much slower relaxation of light particles 
toward equilibrium. A good illustration of this process 
can be found in the top-left panel in Fig. 7. This two- 
scales relaxation mechanism for fi — ► explains the limit- 
ing shape described above for N'(x), see Fig. 1, as well as 
the the positive excess density 5N'(x) in Fig. 2. In addi- 
tion, it also explains the step-like limiting profile for the 
flow velocity: u(x) is the weighted average of the species 
flow velocities Ui(x), the weigh being proportional to the 
species mass, so for /1 — » u(x) coincides with 7*2(2;), 
which is controlled by the fastest relaxation scale. 

It seems also interesting to characterize the region 
where the fluid is far from equilibrium. A canonical 
measure of the scale or typical size associated to this 
region is given by the shock width uj. The shock width 
is usually defined as the inverse of the maximum density 
profile derivative, w(/i) = (^p)max-Q ^ can be also de- 
fined as the standard deviation of the (peaked) distribu- 
tion AN'(x)[l — N'(x)), with A a suitable normalization 
constant. Both definitions give equivalent results. The 
inset to Fig. 8 shows our results for lo(^l). The shock 
width decreases as the mass ratio decreases, scaling as 
w(/i) ~ /1 1 / 3 for fi — > 0.[2a] Therefore nonequilibrium 
effects are more localized around x — for smaller fi. 
However, it is remarkable that despite the shock wave 
gets steeper (uj decreases) as \i decreases, the typical 
scale characterizing relaxation toward equilibrium, A2, 
increases with /i. Hence, although the strong nonequi- 
librium region is reduced as fx — ► 0, the consequent fluid 
evolution to equilibrium slows down. 



IV. CONCLUSIONS 

In this paper we have studied the response of a one- 
dimensional binary fluid to a strong shock wave propa- 
gating into it, on the basis of Boltzmann equation. This 
excitation drives the system far from equilibrium, there- 
fore allowing us to investigate the structure of the fluid 
under nonequilibrium conditions and how equilibrium is 
regained asymptotically. 

Extending to fluid mixtures an elegant method by H. 
Grad, [3113] we have obtained the shock hydrodynamic 
profiles characterizing the transition between the two dif- 
ferent asymptotic equilibrium states. In particular, we 
determine the flow fields as a function of the mixture 
mass ratio /i = m-y/m^ € (0, 1). We find in general that 
all profiles sharpen as (i — > 0. The particle number den- 
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sity field converges in this limit to an asymptotic shape 
characterized by a sudden step-like increase followed by 
a much slower relaxation to equilibrium. The flow ve- 
locity profile converges in turn to a step-like function, 
while the temperature field exhibits a characteristic over- 
shoot which increases with decreasing fi. In addition, the 
density (velocity) of heavy particles behind the shock is 
larger (smaller) than that of light particles, the differ- 
ences increasing as /i — > 0. On the other hand, heavy 
particles are much cooler than light particles right be- 
hind the shock, while the reverse situation holds in front 
of the shock, where heavy particles are now slightly more 
energetic than light ones. 

In order to understand these results, we have per- 
formed an asymptotic analysis of our equations. This 
reveals the emergence of two very different typical length 
scales controlling relaxation as /i — > 0. In this limit, 
the fluid evolves toward equilibrium in twos steps: (i) 
First, heavy particles quickly relax toward their asymp- 
totic downstream equilibrium state, while light particles 
remain close to the upstream state; (ii) at a much slower 
pace, light particles come to equilibrium. In this way, 
light particles ultimately control the fluid's global equi- 
libration. Far behind the shock layer, the hydrodynamic 
fields relax exponentially. Here the multi-scale relaxation 
mechanism is also apparent, with two typical exponen- 
tial scales, Aj, i = 1,2. In particular, Ai, associated to 
heavy particles, goes to zero as /i 1 ^ 2 , while A2, which 
controls the global relaxation to equilibrium, increases 
as fi — > 0, converging to a constant. In addition, the 
size of the strong nonequilibrium region associated to 
the shock layer, as measured by the shock width w(/ii), 
decreases as u)((x) ~ /1 1 / 3 for (1 — > 0. Therefore, the 
fluid's strong nonequilibrium behavior gets more local- 
ized around x — as \i decreases, despite its asymptotic 
relaxation to equilibrium, as given by A2, slows down 
with decreasing \i. 

These results are specially interesting at the light of re- 
cent numerical studies of similar Id binary fluids. 0, 0, 
H5II18I [2(il.l37| In particular, simulations show that when- 
ever the Id binary fluid is perturbed away from equilib- 
rium, the consequent relaxation happens in such a way 
that light particles always tend to absorb more energy 
than heavy ones, as predicted here. In addition, light 
particles are observed to equilibrate more slowly, there- 
fore controlling the system global relaxation. This sup- 
ports the presence of two different relaxation scales for 
/x — ► 0, validating our theoretical results. 

The success of our Boltzmann equation approach to de- 
scribe some of the structural and relaxational anomalies 
observed in the Id diatomic fluid would suggest extending 
the present approach to understand its anomalous trans- 
port properties (see Section I)[I1E1E1E1ESE1E3- 
Unfortunately, this is not possible. Anomalous transport 
in Id is a direct consequence of the long-time power- 
law tails characterizing the decay of correlation func- 



tions in the fluid. The microscopic origin of these tails 
is associated to dynamically correlated sequences of col- 
lisions among fluid molecules, which are not described 
by Boltzmann equation. |38j That is, the main hypothe- 
sis on which Boltzmann equation is based, namely the 
molecular chaos hypothesis. [2ll E2I breaks correla- 
tions present in the fluid which are responsible of the 
power-law tails. 

The strong correlations emerging in Id fluids, and its 
role in the shock wave problem, are interesting issue 
which deserve attention. |37t l3flj | Further work in this di- 
rection, taking into account these correlations in general- 
ized kinetic theories, would be highly desirable to better 
understand the nonequilibrium behavior of simple one- 
dimensional systems. 
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